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By Geobge J. HraoiNs 
SUMMARY 

This paper describes and develops methods by which the aerodynamic characteristics of an airfoil 
may le calculated with sufficient accuracy for use in airplane design. These methods for prediction 
are based on the present aerodynamic theory and on empirical formulas derived from data obtained 
in the N. A. 0. A. variable density wind tunnel at a Reynolds Number corresponding approximately 
to full scale. 

INTRODUCTION 

Since the time that Eiffel first developed the use of the wind tunnel for obtaining the 
aerodynamic characteristics of airfoils, a great amount of data has been accumulated. Many- 
attempts have been made in the analysis of these data to generalize and thereby develop a 
method for the prediction of airfoil characteristics. These have been more or less successful. 
In all cases though, the experimenters have been handicapped by "scale effect." This unknown 
factor undoubtedly contributed a great deal to the differences found between the models and 
the full scale wings. 

In theoretical aerodynamics, the air is treated as a nonviscous fluid. Consequently, the 
profile drag of an airfoil can not be accounted for. However, by assuming that the fluid out- 
side of the "boundary layer" is nonviscous and that within the "layer" is viscous, one can 
combine the pure theory fairly well with the experimental modifications and obtain satisfactory 
results. Scale effect, also a factor caused by viscosity, can not be determined by theoretical 
treatment. 

Empirical derivations from test data modified in accordance with the above theories lead 
to certain formulas which may be used for determining airfoil characteristics. Heretofore such 
expressions have been rather unreliable because of the scale effect mentioned above. Now, 
however, there are available results of tests made under conditions of full scale on airfoils in 
the variable density wind tunnel of the National Advisory Committee for Aeronautics. These 
data then can be used for developing more accurate expressions for prediction. 

This paper covers the development of such expressions from experiments in the variable 
density tunnel and from the present aerodynamic theory. A comparison is made showing the 
differences between predicted airfoil characteristics and observed data from tests in the wind 
tunnel at 20 atmospheres density at a Reynolds Number corresponding to full scale. 

The author acknowledges the suggestions and aid given in the development of this paper 
by E. N. Jacobs and M. Knight, of the Langley Memorial Aeronautical Laboratory. 

THE PREDICTION OP AIRFOIL CHARACTERISTICS 

The main characteristics of an airfoil section are the lift coefficient, profile drag coefficient, 
and the pitching moment coefficient for any angle of attack within the flying range. When 
these characteristics are known, any other derived characteristics may be determined and suf- 
ficient data is then available for use in airplane design. Heretofore, it has been necessary to 
use wind-tunnel tests for the particular airfoil to obtain these characteristics. This has been 
particularly true in regard to the profile drag coefficient because of the lack of theory to cover 
its computation. Test data from the ordinary wind tunnel have been also subject to correction 
for scale effect; but there are now available considerable data from tests in the variable density 
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wind tunnel made at a Reynolds $Tumber equivalent to full scale. Where the term "full scale" 
is used throughout this report a Reynolds Number of approximately 3,400,000 (a VI of 534, or 
a 7#-foot chord at 50 M. P. H.) is meant. 

It is possible to determine the lift and moment coefficients by theoretical computation 
and the profile drag coefficient- by an empirical derivation for any angle of attack up to the 
burble point. The values of total. 0 Df LfB, C>, etc., can then be computed. 

The useful angular range in lift from zero lift to the burble point has been found from tests 
in the variable density tunnel to vary from 18° for thin low cambered sections to 24° for thick 
high-cambered sections. As a normal airplane is limited to aboutrl5° in its range of angle of 
attack from the high-speed condition to that of landing, it may be seen that only in special 
designs does an airplane land or take off at or very near the burble point. Hence,, if an airfoil 
section of normal shape is to be used, the burble point is of secondary importance. 

THE LIFT COEFFICIENT 

The slope of C L curve. — Prandtl, Munk, and others, in their development of the theory of 
lift for an airfoil of infinite span, find the following: 

C^=?27rsin a (1) 

or 

Cx= , 2ira (radians), approx. (Reference 1). 

and 

= _ « 

For wings of finite span having an elliptical span loading and an aspect ratio of g» 

OL*=a< + ot t (3) 

where a fl =absoluta angle of attack measured from the position of zero lift. 

of # = effective angle of attack, i. e., the angle of attack at which an airfoil of infinite span 
would give the same lift coefficient as the airfoil of finite span under consideration. 
From equation (1) above, 

■ a<=induced angle of attack. 
Munk gives for the induced angle of attack 



Substituting for a 4 and a t 
Solving, 



(aa in radians) (6) 

= 0.1Q965 — ^ (a. in degrees) (7) 
Because of the general inefficiency of the wing and because the air is not a nonviscous 



fluid, experiments show that ^~ <2ir 



or 



From this 
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The value of it as determined from many tests conducted at a high Reynolds Number in the 
variable-density wind tunnel is approximately 0.875. 

The induced angle of attack is affected by the shape of the plan form. Glauert gives the 
following expression for the induced angle of attack 



ai-^j? (1 + 0 (Reference 2) 



(10) 



where r is a factor for dete rmining the additional angle caused by the change in span loading 
from that of the elliptical due to the shape of the plan form. For an elliptical plan form r is 
zero. Values for rectangular plan forms of different aspect ratios are given in Table I. 

Table I 

Correction factor* for rectangular airfoils 



1 Aspect 
j ratio, g 


r 




3 


a li 


0. 022 


4 


.14 


.033 


5 


. 16 


.044 


6 


. 18 


.054 


7 


.20 


.064 


8 


.22 


.074 


9 


.23 


.083 



Substituting these expressions of a« and cc { in equation (3) and solving for C Lj one obtains 

C L |£f^ («. iii radians) 

l + lSrCl + r) 



(11) 



or 



C L =- — ° — (*- degrees) 

l+^p— (1 + t) 

The slope of the curve of 0 L plotted against a then becomes 

dO L __ 0.0960 
Ha , , 1.75S 



(12) 



l+^jFCl + r) 



(13) 



For an elliptical wing, or a wing having elliptical span loading, of aspect ratio 6.00, where r is 
zero, 

^=0.0743, 

and for a rectangular wing of the same aspect ratio 

^=0.0714. 

Angle of zero lift. — 

Munk has developed (Reference 3) a method by which the angle of zero lift a La may be 
found. He gives 

+F m | (in degrees) (14) 

where 

<*lo — angle of attack where the lift is zero, measured from the chord line. 
10*397—30 11 
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F a a factor. 

£ = ordinates of the mean camber line at a point* (Z) on the chord line minus the ordinate 

of the mean camber line at the trailing edge. 
c = chord of the airfoil. . . . 
For ordinary use, the 5-point method given is sufficiently accurate. This gives 

-a Lo = 1,252.24 1 + 109.05 ^-+-32.596 |+15.684 ^+5.978 | 5 » in degrees (15) 
c c c c c 

where etc., are the oxdrnates of the mean camber line at the points 

Zi = 99.458% c 

Z f =87.426% c. 

X z = 50.000%" c 

X 4 =12.574% c 

X 6 =0.542%c, 
respectively. • 

By substituting the correct values in the above equation, the angle of zero lift can be 
evaluated. Then, knowing that the lift varies very nearly as a linear function of the angle 
of attack and after determining the slope from equation (13), the curve of 0 L vs. a may be 
readily dr$wn for the section under consideration. u {gee fig. 1.) As stated before, the burble 
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Figure I.—Cl vs. o, hl. must be determined for each airfoil 

point is not known, but this is only of secondary importance. . However, it is of interest to 
note that the highest maximum C L ever recorded, in the variable density wind tunnel on a normal 
airfoil, i. e., without flaps, slots, etc., atra Reynolds .Number equivalent to full scale, is 1.50. 



THE DRAG COEFFICIENT 



The total drag of a wing may be divided into two parts, the profile drag and the induced 
drag, or 

C D =C D ,+ C Di (16) 
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Theory gives 



^i"~jr ( for elKptically loaded wings, Reference 1) 

Ci _ n , ClS 



then 



(17) 



(18) 



It has been usually assumed that G D(f is constant for all values of G L ; but from tests on many 
airfoils at a high Reynolds Number in the variable density tunnel, it appears that this is not 
the case. Mr. Knight, of the laboratory staff, pointed out that the variation in. profile drag for 
these airfoils was similar and seemed to follow a power law- Let 

C D9 +=G DL AC D , (19) 

where Cd Lo is the profile drag coefficient when the lift is zero and when there are no disturbances 
or burbling effects from the lower surface (fig. 2), and whore AGj> 0 is the additional profile drag 
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Ficuri 2 — C;> # vs. Cl 

coefficient, increasing as the lift increases. 

For airfoils that are not extreme in shape and design, values of &> Lo may be found from the 

charts in Figures 3, 4, and 5. These curves of G Blt0 against thickness and against camber 
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Camber, percent chord 

Fig use 3.— Variation of Cn Lm wfta camber 
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Figce£ 4.— Variation of Cd l with thickness 
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Camber, percent chord 
Figure 5.— Additional ACz> due to the effect of camber 

were plotted from data obtained in the variable density tunnel. The accuracy of these curves 
can be surmised by the scattering of points plotted; the error, however, is only a small per cent 
of the total wing drag and still a smaller per cent of the. total airplane drag. 

From Figure 6, where AC Do has been plotted against C L from 22 tests on different airfoils 
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Figure 6.— aCd 0 vs. Cl 



at high Reynolds Numbers, one finds in support of Mr. Knight's suggestion that, as a simple 
approximation 

AC» 0 = 0.0093 C L ™ 5 (20) 
This holds true approximately for all normal-shaped sections. Hence, 



and 



Ql'8 



Total C„ = £^ + 0.0093 C L * M +^p- (elliptical wings) 



(21) 
(22). 



C L *S 



Total C D =C DLo + 0.0093 Cl** + ^t 0- + *) (rectangular wings, Reference 2) (23) 

where a is the factor for the additional induced drag caused by the change in span loading from 
the elliptical. (See Table I.) - 

From the above expressions the polar curves can now be determined up to the burble point. 
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Figure 7 is a chart showing the relationship of the different drag components of a wing or 
airfoil and of an airplane. The different curves in the upper half of the chart show the effect 
of aspect ratio and represent the part of the drag that varies with a change, in the lift; this 
portion of the drag is independent of the choice of the section and is equal to the sum of the 
induced drag and the additional profile drag. The lower part shows the drag dependent on 
the airfoil section and the other parts of the airplane. This chart is convenient for use for 
determining rapidly the characteristics of any airfoil section or any airplane. 




-6 W 1.2 L4 U6 LB 
FiiiURi 7.— Airplane characteristics showing the relation at the different draj components 
LIFT-DRAG RATIO 

Values of L/D can be obtained as usual by dividing the C L by C D and may be plotted 
against the angle of attack or the lift coefficient as desired. Values of a y of course, are deter- 
mined from the curve of C L vs. a. 

THE MOMENT COEFFICIENT 

Reference 3 also gives a method for evaluating the angle of attack when the pitching moment 
about the 50 per cent chord point is zero, a^ 0 



where 



oar, - 62.634 (^-^)» in degrees 



(24) 



£= ordinate of the mean camber line at a poinfr (X) on the chord minus the ordinate of 
the mean camber line at the tr ailing edge. 
X A -95.74% c. 
X fl =4.26%c. 

When the airfoil is in such a position that the moment about the 50 per cent point of the 
chord is zero, the resultant force will obviously pass through this point. Neglecting the moment 
due to the drag force which is very small, the moment about any other point on the chord can 
be computed by obtaining the product of the lift force and its lever arm Q) about the point. 
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(See fig. 8.) By this method, the moment about a point at 25 per cent of the chord. is de- 
termined. Munk shows theoretically (Reference 1) that the moment about this point is con- 
stant for all angles of attack and values of lift. The^moment about the quarter chord point 
becomes 

M&*=-Lfa u ^ Xi=ix~ cos «jjf 0 =Lx ^approximately (25) 
and the moment coefficient is 

M L r ^("u ) G H*m \ 1 

G "4=^-&$ x i~-± or &4i-ry^"i x0 - 0 ^" 0 - 0W - (26) 

where a' is the angle .o^,— «jr § . 0^ is negative giving moment) when Ck(a Vfl ) is positive. 

The C Meh curve may be then plotted by drawing a straight-line parallel to the 0 L axis and with 
a Gx c n value as found above. 

Direction of motion for zero mamgnf fobouf °^°^ } 
f ^-Direction of motion + _ ... 

h<4 




^ Mc/4 " cf l 4 f c/2 

Direction of motion ot zero angle of attack 



FlOUBE 8 

CENTER OF PRESSURE COEFFICIENT 

The center of pressure coefficient can be readily determined from the other characteristics 
already computed by the expression . 

^ = Q:25 - ^cos^ S in« or 0,-0.25-%- , appro*. (27) 

The latter expression is sufficiently accurate for all general purposes. The graph of O v vs. C L 
can be plotted in the usual manner. It is sometime convenient to know the C y in per cent 
chord aft of the quarter chord point. For this 

Q 

0. F.= - 100 -g 1 , per cent chord. (28) 

COMPARISON OF PREDICTED AND OBSERVED CHARACTERISTICS 

Figures 9, 10, and 11 show the characteristics of ap N. A. C, A, M6 airfoil, a Clark Y airfoil, 
and an R. A. F. 15 airfoil, all of aspect ratio 6.85 and with a rectangular plan form. There 
are given on these charts graphs of C Li 0 D) L/D, and Cj^'vs. a, and C D and C p vs. C L . The 
calculated curves are shown as solid lines .and the observed data are shown by points and dashed 
lines. These latter data are from tests at 20 atmospheres density in the variable density wind 
tunnel, a Reynolds Number corresponding to full scale. ■ It is interesting to note the close 
agreement that, is obtained between. the computed .and the observed curves for each of the 
different characteristics.. .Tl^e greatest difference occur? with'the Clark Y section when the 
characteristics are plotted against the angle of attack. This seems to be mainly as the angle 
approaches the critical burbling condition. Since the is computed from the C L) it also is 
at variance with the observed values in this region.^ These discrepancies seem to offset each 
other in the polar curve for there the agreement is good. The predicted L/D curve is slightly 
low mainly because the observed drag of the Clark J is low for a section of its thickness. 

Figures 12, 13, 14, 15, 16, and 17 show the profile drag coefficients for six common sections, 
the U. S. A. 35A, N. A. C. A. M6, U. S. A. 27, R. A. P. 15, Clark Y, and the IT. S. N. P. S. 4, 
plotted against Cz,. Predicted and observed values, are shown here in the same manner. In 
these charts the 0 Do scale is twice the usual one to show the differences more clearly. 
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Fig ubs 8.— Characteristics of M6 airfoil 
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Figure 10 —Characteristics of Clark Y afrfofl 
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Figure 11.— Characteristics of R. A. F. 15 airfoil 
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Figures 12, 13, 14, 15, 16, 17.— Profile drag characteristics 
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CONCLUSION 

The aerodynamic characteristics of an airfoil can be predicted by the above methods with 
sufficient accuracy for use in airplane design. From comparison with observed wind tunnel 
test data at "full scale" it is found that very close agreement is obtained between the predicted 
and the actual characteristics from tests. 



Langley Memorial Aeronautical Laboratory, 

National Advisory Committee tor Aeronautics, 
March U, 1928. 
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TABLE OF SYMBOLS 



C l Absolute lift coefficient, L/qS. 

Ol^ u j Lift coefficient at the angle aji^ 

C d Absolute drag coefficient, D[qS* 

Cd { Induced drag coefficient. 

C D , Profile drag coefficient. 

C d l Profile drag coefficient when the lift is zero. 

AC o m Additional profile drag coefficient or € D m =C d Lo + AC d 0 . 

C Dp Parasite drag coefficient. 

Cu efi Pitching moment coefficient, -^§f* 

C p Center of pressure coefficient as a fraction of the chord. 

C. P. Center of pressure. 

L Lift. 

D Drag. 

M cfi Pitching moment about the quarter chord point. 

L/D Lift-drag ratio. 

a K Absolute angle of attack, measured from the position where the lift is zero. (See fig. 1.) 

cti Induced angle of attack. 

a. Effective angle of attack. 

a Lm Angle of attack when the lift is zero, measured from the chord line. 

ajt m Effective angle of attack when the moment about the half chord point is zero. 

a* Angle of attack equal to a^— ajf tt . 

b Span of airfoE. 

c Chord of airfoil. 

S Plan form area of airfoil. 

q Dynamic pressure, }pF*. 

P Density of air. 

V Velocity. 

ifc Correction factor for inefficiency of airfoiL (See text.) 

t Correction for additional induced angle of attack for rectangular wing. 

cr Correction for additional induced drag coefficient for rectangular wing. 

ht etc. Ordinates of mean camber line of airfoil at points (X) on chord; for use in deterniining a i# . 

$U» Zb Ordinates of mean camber line of airfoil at points (X) on chord; for use in determining cejf,. 
X tf Xt, etc. Points on the airfoil chord. (See text.) ■ 

Xa, X B Points on the airfoil chord. (See text.) 
fu fa ete. Multiplying factors. (See text.) 



